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ABSTRACT 



We have developed a data reduction procedure to extract multiple spectra from a single two- 
dimensional Space Telescope Imaging Spectrograph (STIS) image of a crowded stellar field. This 
paper provides a description of our new technique, utilizing a STIS ultraviolet spectral image, 
acquired with the G140L grating and the 52" x 2" aperture, sampling a concentration of O 
■ and B stars in the central region of the NGC 604 starburst in M33. The software routines can 

disentangle and produce reliable ultraviolet spectra of stars with angular separations as small as 
0.055". Use of the extraction slit, based upon our model of the spectral cross dispersion profile, 
generates spectra with slightly higher resolution than the STScI standard processing. Our results 
clearly show that the spectral imaging capability of STIS represents a powerful tool for studying 
' luminous stars in the star-forming regions of the Local Group. 



Subject headings: ISM: individual (NGC 604) — instrumentation: spectrographs — techniques: 
image processing — instrumentation: high angular resolution — techniques: spectroscopic 



1. INTRODUCTION 

One of the designed functions of the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble 
Space Telescope (HST) is to obtain two-dimensional stellar spectroscopy similar to objective prism spec- 
troscopy from the ground. However, the high angular resolution capabilities of the HST and STIS are far 
better than can be obtained by comparable instrumentation from ground-based observatories. Theoretically, 
spectra of stars with separations of less than 0.1" can be disentangled and reliably studied. The STIS has 
another advantage in that it can perform ultraviolet through optical (1150-10,000 A) spectroscopy of indi- 
vidual stars in star-forming regions anywhere in the Local Group. Despite these instrumental capabilities, no 
significant two-dimensional studies have been made of crowded stellar fields using the STIS. This has been 
in part because of the lack of any data reduction and analysis software capable of extricating the overlapping 
spectra of stars in a densely populated field-of-view of a STIS spectral image. 

In this paper, we describe our data reduction software developed for processing STIS spectral images 
of crowded stellar fields acquired with the G140L grating and the large 2" wide aperture. Specifically, 
we present a description of the software and techniques used in reducing and analyzing a STIS ultraviolet 
spectral image of the OB association, NGC 604. 
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We observed with STIS the UV-bright concentration of and B stars in NGC 604, the luminous star- 
forming region located in the spiral arms of the nearby galaxy M33. NGC 604 is one of the brightest H II 
regions in the local Universe. Furthermore, its close proximity (d= 840 kpc) makes it an excellent target to 
probe the details of enhanced star formation activity. Our objective was to observe and resolve the young, 
embedded stars in the luminous H II region and obtain ultraviolet spectra of large numbers of O and B 
stars from a single image. A full discussion of the subsequent analysis of these data are presented in the 
companion paper, Bruhwcilcr, Miskey, & Smith Neubig (2003; Hereafter Paper 2). 

When obtaining an image with multiple objects in a single exposure, the overlapping spectra provide 
particular data reduction challenges. This is certainly the case for our STIS data of NGC 604. In what 
follows, we describe aspects of data acquisition that affect the data reduction and discuss the basic technique 
used to produce accurate flux and wavelength calibrated spectra. We further present examples that give an 
assessment of the data quality of the extracted spectra. 

2. OBSERVATIONS 

We obtained a two-dimensional, ultraviolet, spectral image (dataset ID: o4xl01040) of the brightest 
concentration of stars in NGC 604 utilizing the STIS aboard HST on 1998 August 27. A 1720 second exposure 
was acquired with the FUV-MAMA detector, G140L grating, and 52" x 2" aperture (that effectively samples 
a 25" x 2" field-of-view (FOV) for the MAMA detector). The aperture was aligned with the position angle 
of the image y-axis at P.A. = -165.009 degrees East of North to maximize the number of bright stars in the 
slit, based on previous pre-COSTAR HST Wide Field Planetary Camera (WFPC) imagery, as presented in 
Drissen, Moffat, & Shara (1993), and Wide Field Planetary Camera 2 (WFPC2) imagery (see Table 1). 
See Figure 1 for an overlay of our STIS aperture on an archival WFPC2 image taken with the F336W filter. 
(Also, see Paper 2.) For stars positioned along the long axis of the aperture, the low resolution (~ 2-3 A) 
G140L grating provides an effective wavelength range from 1150 to 1736 A. The plate scale of the 1024 x 
1024 pixel image is 0.0244 arcsec/pixel for the FUV-MAMA. 

Although spectral images can be obtained of a larger number of stars using the slitless mode, the 
geocoronal emission would contaminate the entire FOV, greatly increasing the background and lowering the 
signal-to-noise. The 2" wide aperture is wide enough to include a large fraction of the O and B stars in 
NGC 604, but still restrict the region of the spectrum affected by the high Lya background. See Section 3.4 
for more discussion of the astronomical background. 

Our data show at least 60 UV spectra of O and B stars (see Figure 2). Geocoronal H I- Lya emission 
outlines the full 25" length of the MAMA FOV of the 2" wide STIS aperture. Each spectrum clearly shows 
sharp interstellar absorption produced by gas in both the Milky Way and M33 galaxies. These components 
are unresolved spectroscopically, since the radial velocities of M33 and NGC 604 are v r = —179 km s _1 and 
—226 km s -1 , respectively (Helou et al. 1995), and appear as a merged single velocity component at the G140L 
resolution. 2 These sharp interstellar lines include Si II AA1190,1193, N I AA1198,1199,1200, Si III A1206, 
Si II A1260, O I A1302, Si II A1304, C II AA1334,1335, Si IV AA1393,1402, Si II A1526, C IV AA1548,1550, 
Fe II A1608, and Al II A1670. The interstellar Si IV and C IV features are superposed upon broad stellar 
P Cygni mass loss features. Other prominent broad features include interstellar H I-Lya, as well as stellar 



2 The radial velocities are from the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the NASA. NED is available at: http: / / nedwww.ipac.caltech.edu. 
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P Cygni N V AA1238,1240 and He II A1640 emission. The actual x-positions on the MAMA detector of the 
spectral lines are shifted, corresponding to the x-positions of the stellar sources in the 2" wide aperture. 

Besides the HST imagery of NGC 604, we incorporated STIS spectral images of other individual stars 
to both design and test our data reduction software. See Table 1 for a list of the relevant datasets. 

3. DATA REDUCTION 

3.1. MODELING THE STIS CROSS DISPERSION PROFILE 

In this section, we describe the process for modeling the cross dispersion profile (CDP) of a STIS 
spectrum. To facilitate our discussion, in Figure 3 we present a flowchart summarizing the steps. Because 
of the many subroutines and data products used in our data reduction, we also provide Appendix A, which 
contains a glossary of our terms. The reader should feel free to refer to Appendix A as necessary. We give 
a more detailed discussion of how and why we perform the CDP modeling below. 

To separate closely spaced spectra on the MAMA detector face of the STIS, it is first necessary to 
accurately model an isolated spectrum produced by the MAMA/G140L configuration. This will generate 
key parameters defining the shape and variation across the detector of the cross dispersion profile of the 
spectrum. We have selected a 1024 pixel x 1024 pixel STIS/G140L spectral image of the bright, nearby, 
hot, white dwarf, GD153 (ID: o4a502060; See Table 1). We choose the UV spectrum of this star to model 
because GD153 represents a STIS calibration target and is a continuum source with few photospheric or 
interstellar lines. 

All STIS/G140L datasets used in this study undergo preliminary processing in advance to produce a flat- 
fielded image. This processing includes statistical error initialization, conversion to count rates, subtraction 
of a dark rate image, division by a flat-field, and bad/hot pixel processing. 

Seven vertical crosscuts, SLICEi(y), are taken of the white dwarf spectrum in the flat-fielded, background- 
subtracted STIS image. The subscript i denotes the index (1-7) corresponding to a specific x-value, and the 
y represents a 1024 element array perpendicular to the dispersion. These seven crosscuts are each the sum 
of seven pixel wide slices centered at Xj=[69, 290, 511, 732, 953, 995, 1017]. These x-positions correspond to 
default, unshifted wavelengths of A;=[ 1164 A, 1293 A, 1422 A, 1551 A, 1680 A, 1705 A, 1718 A]. The spac- 
ing between these vertical slices at longer wavelengths is smaller because the response curve of the MAMA 
detector changes more rapidly at higher x-values (longer wavelengths). (See Figure 4.) Each of these slices 
is fit very accurately with a function represented by the sum of three Gaussians, YFITi(y), each centered 
on the peak of the spectrum. See Figure 5 for three of the CDP crosscuts (SLICE 2 , SLICE 3 , SLICE4) 
overplotted with heavy dotted curves denoting our triple-Gaussian fits (YFIT2, YFIT3, YFIT4). Each 
YFITi(y) is then normalized by dividing by the total integrated area of its triple-Gaussian profile. These 
seven normalized fits, NORMi(y), determine the shape of the stellar spectrum perpendicular to the disper- 
sion (CDP) as a function of x (wavelength). These normalized triple-Gaussian profiles represent our new 
weighted extraction slit and replace the box-car slit used in standard STIS data processing. 
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3.2. DEVELOPING OUR DATA REDUCTION PROCEDURE 

In this section, we describe the development of our data reduction procedure, STISFIT, written by C. L. 
Miskey in the Interactive Data Language (IDL). Wc utilize this set of programs for the subsequent spectral 
extractions. The simple flowchart in Figure 6 outlines the development process for which we provide a more 
detailed discussion below. 

We use our software first to perform the background subtraction for the region in the STIS image 
surrounding the stellar spectrum of GD153, which is trivial for data with a single star in the aperture. Two 
background slices (one above and one below), each 5 pixels wide by 1024 pixels long, are taken from areas 
immediately adjacent to the spectrum, but far enough removed not to be contaminated by stellar flux. This 
background array is averaged and then subtracted from the data. 

Second, a generic template of a STIS spectral image is constructed using the NORMi(y) arrays. A 
two-dimensional template array is built, FIT(x,y), 1024 pixel columns by 1024 pixel rows. Each i of 
the NORMi(y) is associated with a proper wavelength. Since GDI 53 is centered in the aperture, we use 
the default, unshiftcd, heliocentric wavelength scale, W{x), spanning 1150-1736 A. Each x (or wavelength) 
column of the template is individually created by linearly interpolating between corresponding values in the 
two bracketing vertical fits, NORM i ^ 1 (y) and NORM^y). For values of x<69 (i=l), NORM 1 (y) is used. 
Likewise, FIT(x > 1017, y) = NORM 7 (y). 

The previous step addresses variations in the CDP across the array in x. However, we must also consider 
corrections in the y-direction. 

A STIS stellar spectrum is not straight, but has a curvature that varies along the dispersion. This 
curvature is not fixed, but also varies relative to y-position on the detector. This spatial distortion has 
been accurately delineated by Don Lindler, who has compiled a reference file, exttab, for the G140L grating 
containing distortion curves for 27 different positions along the y-axis (Lindler 1999). Similar exttab files 
assembled by Lindler exist for all of the other STIS gratings. The curvatures for y-positions intermediate 
to those positions in the table are calculated through linear interpolation. We define the y-position of a 
spectrum, represented by the variable cent in our software, as the y-location of the center of the spectrum at 
x=512 in the 1024 x 1024 pixel STIS spectral image. This definition is used throughout the paper. Using the 
G140L exttab file and the y-location of GD153 on the detector, the template array, FIT(x,y), can easily be 
corrected for spatial distortion. The location of the spectral peak in FIT(x, y) is mapped to the calculated 
distortion for each column of the array. The spatial distortion across the STIS detector for three different 
y-positions is plotted in Figure 7. A horizontal, dotted line runs through the star's central location at x=512. 
Note that the deviations in y of the spectrum are most pronounced at extreme low and high values of x. 

Finally, this corrected data template is applied to the STIS image for GD153. Each of the 1024 columns 
of the FIT(x,y) array is fit to each of the 1024 columns of the flat-fielded, background-subtracted data, 
IMAGE(x,y), using a linearization of the fitting function, to produce an array containing the model spec- 
trum, FINALFIT(x, y). The routine generating FINALFIT(x, y) uses the same technique as Bevington's 
CURFIT (Bcvington 1969). Using the notation of Bevington, for a function Y (the column of FIT(x,y) 
being fit) of parameters A(I), the routine uses a linear expansion of Y as a function of A(I). The array A 
has n values (where n equals the number of stellar sources; here n=l) denoting the amplitude(s) of the 
triple- Gaussian function(s). The Y is evaluated for a set of parameters A and the array of independent vari- 
ables X. The procedure uses generalized least-squares to find the parameter increments which minimize x 2 . 
The program automatically iterates through parameter step sizes, which are controlled by input parameters. 
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Standard deviation errors are approximated by the inverse of the curvature matrix, (l/cr 2 ). 

The resulting fits show extremely good correlation. We calculate the least squares fit between the actual, 
background-subtracted image and our STISFIT model image. 

T,l=f T,T=- 50 [IMAGE(x, cent + y) - FINALFIT(x, cent + y)} 2 

Least squares = 77^5-, ^ x 100%. 

E^i -bo [IMAGE(x, cent + y)] 2 

For GD153, a result of 0.39% is obtained. 

The STISFIT net spectrum, SPECTRA(x), is then created by summing an eleven pixel wide hori- 
zontal slice of the FINALFIT(x,y), centered on the stellar y-location. This extraction is not a straight 
horizontal cut, but follows the curvature on the detector characterized by the exttab file, as described above. 
SPECTRA(x) is in units of counts second -1 pixel -1 . 

We derive a new sensitivity curve, S(x), by requiring that our spectral extraction yield the same mea- 
sured fluxes versus wavelength for GD153 as the CALSTIS extraction. CALSTIS is the standard STIS 
processing routine used in the Laboratory for Astronomy and Solar Physics at NASA/Goddard Space Flight 
Center (Lindler 1999). Also, a very similar version of CALSTIS is used in the STScI pipeline processing 
(Brown 2002). Our STISFIT sensitivity curve, S(x), is determined by dividing the GD153, flux-calibrated 
spectrum generated by CALSTIS by our net STISFIT spectrum, SPECTRA(x). This curve is dependent 
on wavelength. 

With x corresponding to the unshifted wavelength array, W(x), the STISFIT sensitivity curve is 
multiplied by the STISFIT generated net spectrum to yield the final, flux-calibrated, stellar spectrum; 
S(x) x SPECTRA{x) = FLUX(x). The flux spectrum is in units of ergs sec" 1 cm" 2 A" 1 . 

Upon inspection of Figure 8, one sees that our final flux spectrum of GDI 53 differs slightly from the 
CALSTIS spectrum at low wavelengths. This is because the MAMA detector has low sensitivity in this 
region, which leads to low signal-to-noise for the data and to a correspondingly less accurate sensitivity 
curve in the region shortward of 1200 A. Therefore, we have chosen to smooth the derived S(x) through this 
region to reduce the effects of noise and produce a more reliable sensitivity curve at these wavelengths. 



3.3. ASSESSMENT OF RESULTANT SPECTRA 

To test how closely our STISFIT procedure reproduces CALSTIS spectra, we calculate the percentage 
difference and average deviation between the extracted flux spectra generated by these two methods. Here, 
the percentage difference and average deviation are given by: 

Et-i 4 [CALSTISix) - FLUX(x)} 

Percentage difference = L ^ x ~ 1 1 ln9A ^ x 100% 

E'li CALSTIS(x) 

and 

, . . Et-i 4 J[CALSTIS(x) - FLUX(x)} 2 ^ nnM 

Average deviation = ^ x ~ 1 v 1 in9/l — x 100%. 

Ex=i CALSTIS(x) 

We use the STIS/GI40L datascts of three stellar objects, GD153, NGC 346-368, and Sand 3. All arc 
stellar point sources, but each has a different type of spectrum. As previously stated, GD153 is a hot, white 
dwarf, continuum source with few photosphcric or interstellar lines. The object NGC 346-368 is a luminous 
O star in the Small Magellanic Cloud with a rich photosphcric absorption spectrum, while Sand 3 is a hot, 
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pre-white dwarf with a Wolf-Rayet emission line spectrum. (See Figure 8.) For the percentage differences, 
we obtain -0.06%, 2.34%, and 1.20%, respectively. The percentage difference indicates how closely the total 
flux is reproduced between the CALST1S and STISFIT spectra. The small values indicate that the total 
flux is conserved to within 2.5% in our worst case. These differences are in accord with the variations in 
sensitivity of the MAMA/G140L (Stys & Walborn 2001) 3 . We expect a very small value for GD153, since 
it is the spectrum used to construct the STISFIT procedure, and one of our requirements of S(x) is that 
STISFIT yield the same flux values as CALSTIS. The insignificant difference is due to smoothing of the 
STISFIT sensitivity curve at lower wavelengths (A < 1200 A). The largest variation is for an object with a 
rich line spectrum. For the average deviation, we calculate 0.64%, 3.77%, and 2.82%, respectively. These 
results give us sufficient confidence in our procedures to apply them to our NGC 604 dataset. 

3.4. APPLYING OUR DEVELOPED MODELS AND SOFTWARE TO THE STIS 

DATASET FOR NGC 604 

We apply the newly developed routines, STISFIT, to our STIS/GI40L spectral image of NGC 604. This 
set of programs has the capability to produce extracted, flux and wavelength calibrated, spectra for clusters 
of stars as well as isolated stars. The final calibrated spectra are produced in a two step process. We include 
a simple flowchart in Figure 9, which outlines the reduction procedure described in detail below. 

The first step is to use STISFIT to create an extracted, uncalibrated spectrum for each star. We 
start by removing the background from the flat-fielded STIS image. The sky background has several major 
contributions. The geocoronal emission is comprised of the diffuse H I-Lya at 1216 A and the much weaker 
O I Af 304. These geocoronal emission lines fill the entire 2" wide slit in the spatial dimension and yield two 
vertical bands, one bright and the other much fainter, in the spectral image. The H I-Lya band is quite 
distinct in the spectral image displayed in Figure 2. The weak O I emission is at most 10% of the H I-Lya 
in brightness, depending on the time of the observation and the position of the target relative to the Earth. 
In our observation, O I is negligible. 

Because the scattered light from dust is essentially scattered UV stellar flux, its spectral characteristics 
mimic the spectra of the adjacent stars. Consequently, care must be taken in applying any background 
subtraction. This is especially the case for stars with low UV flux, where the accuracy of the background 
subtraction has a larger effect on the net spectrum. In the UV wavelength range sampled by the G140L, 
nebular emission lines appear to be a negligible contribution to the background. 

The most variable procedure in our data reduction scheme is the background subtraction. Ideally, one 
wants to determine the background as close to the stellar source as possible. Depending upon the degree of 
crowding in the FOV, different techniques may need to be applied. As long as there is a reasonable amount 
of area free of stellar flux in the image, the standard background subtraction, averaging two bordering back- 
ground strips, can be performed. However, for a spectral image that contains very little area uncontaminated 
by stellar flux, a background image must be tailor-made. Our NGC 604 image presents such a severe case. 

In our NGC 604 spectral image, the densely populated stellar field makes it impossible to obtain direct 
background measurements adjacent to most of the individual spectra and subtract the background in the 
usual manner. Hence, we create a 1024 x 1024 pixel array of the underlying background to be subtracted. To 



3 STIS Instrument Science Reports are available at: http://www.stsci.edu/hst/stis/documents/isrs. 
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map the background relative to wavelength, two horizontal strips of the STIS image (one in the lower portion 
and one in the upper portion of the dataset), 1024 pixels long by 10 pixels wide, are extracted and averaged 
from regions of the image containing no detectable UV stellar sources. The section containing the geocoronal 
H I-Lya emission is assigned zero weight, and a 6th order polynomial is fit along x to the averaged background 
strip. The H I-Lya band is then modeled with 2nd order polynomials for the regions corresponding to the 
edges of the 2" wide aperture and a mean constant value denoting a flat plateau for the top. These curves are 
combined to form the model, continuous, functional fit of the background in the x-direction, BX{x). Several 
vertical crosscuts of the image are then taken, and the minimum background values in the y-direction arc 
measured to create an array of background levels perpendicular to the dispersion, BY(y). The two arrays 
BX(x) and BY(y) are combined, yielding a complete two-dimensional model, B(x,y), for the background 
of the STIS spectral image of NGC 604. We find this method of background determination quite reliable. 

After the background model is subtracted from the image, STISFIT begins the spectral extractions. 
It is very important to accurately determine the location of the center of each stellar spectrum, cent(n), 
preferably on the sub-pixel level. This is necessary to properly apply our triple-Gaussian spectral extraction 
algorithm. Isolated spectra are extracted individually for faster processing. For regions where a group of 
stellar spectra overlap or contaminate the local background, these spectra are all modeled together. The 
variable n equals the number of stars in each extraction set. For each grouping, a FIT(x,y,n) array 
is assembled with n triple-Gaussian functions representing each of the point sources, corrected for stellar 
locations in the aperture and spatial distortions on the detector face. Next, the template FIT(x,y,n) is fit 
to the flat-fielded, background-subtracted IMAGE(x,y). A model for each spectrum, SPECTRA(x,n) 7 is 
then extracted from the grouped FINALFIT(x,y,n). 

The second step is to calculate the wavelength offset of each star and create the final, extracted, wave- 
length and flux calibrated, spectra. First, strong interstellar lines are identified from the uncalibrated spectra. 
These sharp, identified, spectral features serve as fiducial reference wavelengths. They are used to calculate 
shft(n), the x-position (wavelength) offset relative to the default, where the star is centered in the aper- 
ture. Our corrected, heliocentric, wavelength array, W\(x,n), is determined by using shft(n) to shift the 
initial wavelength array, W(x), for each spectrum accordingly. The S(x) array is mapped onto the new 
wavelength scales, producing the corrected sensitivity curves, S\(x,n). The A subscript denotes that the x 
array correlates to a shifted wavelength scale. 

The shft(n) array is then used in a second run of our fitting routine, STISFIT, to properly align 
the NORMi(y) arrays with their corresponding shifted wavelengths when the initial FIT(x,y,n) tem- 
plate is built. The shifted FIT(x,y,n) is fit to IMAGE(x, y), yielding FIN ALFIT(x, y, n). The new 
SPECTRA(x, n) is extracted. The shifted sensitivity curves, S\(x, n), are multiplied by the results from the 
second pass of STISFIT, SPECTRA(x, n), yielding the final, wavelength (A) and flux (ergs sec" 1 cm~ 2 A -1 ) 
calibrated, stellar spectra, FLUX(x,n). 

4. COMPARISONS OF OUR DEVELOPED PROCEDURE WITH STANDARD STIS 

PROCESSING 

4.1. THE SPECTRAL RESOLUTION OF EXTRACTED SPECTRA 

Upon close examination of Figure 8, one finds that our STISFIT extraction, utilizing the triple-Gaussian 
fit to the cross dispersion profile as represented in Figure 5, produces deeper, sharper, absorption features 
than the CALSTIS extraction. These differences are most apparent in the comparison of the processing 
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of the G140L data for the bright O star, NGC346-368, in the H II region, NGC 346, of the SMC. This is 
displayed in the middle panel of Figure 8. Likewise, narrower, sharper N V A1240, Si IV A1400, C IV A1550, 
and He II A1640 profiles are also seen in the STISFIT spectrum in the comparison for the strong emission 
line object, Sand 3, in the bottom panel of Figure 8. 

The slightly sharper features in our STISFIT extractions imply that STISFIT yields improved spectral 
resolution over the CALSTIS processing. This might be expected, since our triple- Gaussian function gives 
higher weight to data near the center of the spectrum, whereas the box-car extraction used in CALSTIS 
gives equal weight to all datapoints. The CDP can be thought of as the total contribution of the individual 
point spread functions (PSFs) produced at every wavelength along the spectrum. When using a narrow 
slit, the PSF contribution is limited to the truncated portion that passes through the slit and is imaged 
on the detector. This restricts the effects of the PSFs at greater distances in the dispersion direction. The 
contribution from non-local PSFs at any point on the CDP away from the spectral peak, perpendicular to 
the dispersion, drops off dramatically due to truncation from the narrow slit. However, the spectral images 
discussed here were obtained with the 2" wide slit. The PSFs are not truncated, unlike the case for the narrow 
slit. Thus, at any specific wavelength, A Q , as one moves away from the center of the spectrum, perpendicular 
to the dispersion, there are higher fractional contributions to the flux from wavelengths farther from A . Since 
our triple- Gaussian extraction slit gives higher weight to points near the center of the spectrum, the derived 
STISFIT fluxes at A are determined by points much more local to A than for the box-car extraction fluxes. 
Other similar diffraction effects are discussed in Chapter 13 of the STIS Instrument Handbook (Lcitherer et 
al. 2001) and in Bowers & Baum (1998). 

To quantitatively determine if STISFIT produces higher resolution data, one can compare both low 
resolution G140L STISFIT and CALSTIS extractions with a higher resolution E140M CALSTIS extraction 
for the same object. The STIS/E140M echelle grating (R ~ 10,000) spans the same wavelength range 
as the G140L. The better G140L extraction should yield spectral line measurements more consistent with 
those obtained from the E140M spectrum. For our comparison, we select the star NGC 346-368, since it 
has been observed with both gratings. We compare the three extractions, and, as expected, the E140M 
data clearly show much sharper spectral features than those of both G140L extractions. We measure the 
equivalent widths of the strong interstellar lines in all three extracted spectra. For the G140L data, the 
STISFIT spectrum yields larger equivalent widths than the CALSTIS spectrum. The STISFIT equivalent 
widths agree better with those measured in the E140M dataset, but are still systematically smaller. This 
comparison indicates that the STISFIT extracted spectra provide better quantitative measurements for 
critical spectral parameters such as equivalent width. 

We also estimate the actual spectral resolution for the two G140L extractions. We do this by convolving 
(smoothing) the E140M data with a single Gaussian profile and correlating the result to the STISFIT and 
CALSTIS datasets. The full width half maximum (FWHM) of the Gaussian smoothing is varied to obtain a 
least-square minimum (minimizing Ecr 2 ). The best matches to the STISFIT and the CALSTIS spectra are 
achieved if the E140M spectrum is convolved with Gaussians with FWHMs equal to 2.2 A (R ~ 640) and 
3.0 A (R ~ 470), respectively. (See Figure 10.) Since the unsmoothed, intrinsic resolution for the E140M 
grating is a factor of 15 to 21 times higher, these two deduced values are good estimates of the spectral 
resolution near 1380 A. Once again, our analysis suggests a higher resolution for the STISFIT extraction. 

All resolution tests were performed on datasets obtained with the 2" wide aperture. We have not used 
our developed software to reduce any narrow slit STIS data. Hence, we do not know if, given the effects of 
PSF truncation, our STISFIT technique would still yield a higher resolution in such cases. 
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4.2. CLOSELY SPACED & OVERLAPPING SPECTRA 

Perhaps the most important feature of our STISFIT process is its ability to extract several closely 
spaced spectra concurrently. By modeling each stellar spectrum with its own triple-Gaussian fit, we are 
able to successfully separate out each star's flux contribution and minimize cross-contamination between 
neighboring spectra. The standard STIS processing method can only extract one spectrum at a time with 
an eleven pixel wide box-car extraction slit. Hence, it has no means of dealing with multiple overlapping 
spectra. 

One excellent example in our NGC 604 dataset is the stellar complex near y=578, composed of three 
overlapping stars, 578A, 578B, and 578C. (See Figure 2.) Their spectral centers, cent(n), are located at 
y=[578.700, 580.945, 583.626] in the image, corresponding to separations of 2.245 and 2.681 pixels, or 0.055" 
and 0.065", in the y-dircction. The STISFIT procedure produces three clean, low-noise spectra with very 
little cross-talk. (See Figure 11.) Moreover, the spectral features and shape are consistent with what is 
expected for 09 stars. This is even the case in the region beyond 1600 A, where the sensitivity of the 
MAMA falls off sharply. 

5. DISCUSSION 

We have presented the basic techniques that compose our STISFIT extraction procedure for STIS 
spectral images of crowded stellar fields. The foundation of our process is a triple-Gaussian extraction 
slit that gives higher weighting to flux values near the center of the spectrum. As illustrated in Figure 
5, the triple-Gaussian profile very closely fits the CDP, which enables us to determine the amount of flux 
corresponding to each stellar source for a group of overlapping spectra. We have made use of the G140L 
spatial distortion reference file such that our extraction procedure can accurately track a spectrum at any 
position in the two-dimensional STIS spectral image. The resulting STISFIT spectra yield deeper stellar 
absorption features. Our tests indicate that the spectral resolution for STISFIT spectra are slightly better 
than that of the standard CALSTIS or STScI pipeline processing. Furthermore, the strengths of the STISFIT 
extraction technique have allowed us to produce reliable, flux and wavelength calibrated, G140L spectra for 
stars separated by as little as 0.055" in the direction perpendicular to dispersion. This is at or near the 
spatial resolution limit of STIS. Additional examples of our extracted spectra can be found in Paper 2. 

We have recently expanded the STISFIT software to include the other low resolution grating configu- 
rations of STIS. Our software routines are now applicable to spectral images obtained by both MAMA and 
CCD detectors, several different gratings, and the 52" x 2" and slitless aperture modes. These improve- 
ments allow us to now obtain continuous spectral coverage, from ultraviolet through optical wavelengths, of 
individual stars in crowded fields. We will present observational results and assessments of these techniques 
in following papers. 

We wish to acknowledge the assistance and invaluable advice provided by Don Lindler, Nicholas Collins, 
Eliot Malamuth, and Charles Bowers; all arc members of the STIS Instrument Development or Support 
Groups at NASA/Goddard Space Flight Center. We also acknowledge support through NASA grant NAG5- 
3378. The NGC 604 STIS spectral image was obtained under Guaranteed Observing Time awarded to Dr. 
Albert Boggess. 
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A. Glossary 

B(x,y) Two-dimensional model array of the background for our STIS spectral image. 

BX(x) Model, continuous, functional fit of the background in the x-direction for our STIS dataset. 

BY(y) Array of background levels perpendicular to the dispersion for our STIS dataset. 

CALSTIS Standard STIS data reduction procedure used at LASP, NASA/GSFC and STScI. 

CALSTIS(x) Flux-calibrated spectrum from a STIS image generated with CALSTIS. 

CDP Cross dispersion profile. 

cent(n) Y-location in pixels of the center of n spectra at x=512 in a 1024x1024 pixel STIS spectral image. 
exttab Reference file containing spatial distortion curves for STIS gratings. 

FINALFIT(x,y,(n)) Resultant array containing the model spectra created by fitting our STIS template, 
FIT(x,y, (n)), to the spectral image, IMAGE(x,y). 

FIT(x,y,(n)) Generic template of a STIS spectral image constructed with NORMi(y), cent(n), and exxtab. 
Accounts for variations in both x and y on the STIS detector. 

FLUX(x,(n)) Final, flux-calibrated, stellar spectra from a STIS image generated with our STISFIT pro- 
cedure in units of ergs sec -1 cm -2 A -1 . 

IMAGE(x,y) Flat-fielded, background-subtracted, STIS image. 

NORMi(y) Seven normalized YFITi(y) arrays, each a triple-Gaussian function delineating the cross dis- 
persion profile of the STIS detector at a specific value of x. Represents our new weighted extraction 
slit. 

S(x) Our default STISFIT sensitivity curve. 

S\(x,n) Shifted STISFIT sensitivity curves, corrected for stellar offset relative to default in the aperture. 

shft(n) The x-position (wavelength) offset relative to default, where the star is centered in the aperture, 
for each n. 

SLICEi(y) Seven vertical crosscuts, perpendicular to the dispersion, extracted from the STIS image of 
GD153. Each is the sum of seven columns and corresponds to a specific x (wavelength) value. 

SPECTRA (x,(n)) Net spectra from a STIS image generated with our STISFIT procedure in units of 
counts second -1 pixel -1 . 

STISFIT Our set of data reduction routines, written in IDL, which extract flux and wavelength calibrated 
spectra from a crowded STIS field-of-view. 

W(x) Default, heliocentric, wavelength scale. For STIS/G140L, it spans 1150-1736 A. 

W\(x,n) Shifted, heliocentric, wavelength array, corrected for stellar offset relative to default in the aper- 
ture. 

YFITi (y) Seven resultant arrays from the fit of triple-Gaussian functions (the sum of three Gaussians, each 
centered on the peak of the spectrum) to the SLICEi{y) arrays. 
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Fig. 1. — Outline of the 2" wide aperture from our STIS observation superposed upon WFPC2/F336W 
imagery of NGC 604. The position angle of the aperture is -165.009 degrees East of North. 
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Fig. 2.— Our STIS/G140L spectral image of bright O and B stars of NGC 604 with stellar complex 578 
denoted. The vertical band to the left in the spectral image delineates the H I-Lya geocoronal emission. 
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SLICEfr): 

Extract 7 vertical crosscuts from STIS image 
ofGD153. (Note: / represents the index 1 -7.) 

* z 

YFITJy): 

Fit each SLICEJy) with a triple-Gaussian function. 

NORMjfy): 
Normalize YFlTfy) arrays. 



NORM! (y) arrays characterize the cross 
dispersion profile of a STIS/G140L spectral image. 



Fig. 3. — Modeling the STIS cross dispersion profile with a MAMA/G140L spectral image of white dwarf 
GD153. (See Section 3.1 for further detail.) 
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Fig. 4. — Seven vertical crosscuts, SLICEi(y), are shown on the STIS/G140L spectral image of GD153. 
The 7-pixel wide by 1024-pixel long slices, centered at Xj=[69, 290, 511, 732, 953, 995, 1017], are summed in 
x. 
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Fig. 5. — Cross dispersion profiles of GD153 overplotted with our triple-Gaussian fits centered at three 
different x-locations. The STIS data are marked with the solid lines. Our triple-Gaussian functions are 
marked with the heavy, dotted lines. 
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IMAGEfcy): 
Flat-field and background-subtract STIS image. 



FIT(x,y): 

Construct 2-D template of STIS image. 



Linearly interpolate between 
NORMJy) arrays to build FIT(x,y). 
Accounts for variations in x. 



Correct FIT(x,y) for spatial 
distortion using cent and exttab. 
Accounts for variations in y. 



FINALFITfcy): 

Fit FIT(x,y) to IMAGE(x,y). 



SPECTRA(x): 

Extract net spectrum 
from FINALFIT(x,y). 



IF 



S(x): 

Derive sensitivity curve. 



FLUX(x): 

Flux -calibrate final spectrum (ergs s _1 cm -2 A 1 ) 
by multiplying SPECTRA(x) by S(x). 



Fig. 6. — Developing our data reduction procedure for a STIS/G140L spectral image utilizing the GD153 
dataset. (See Section 3.2 for further detail.) 
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Fig. 7. — The spatial distortion of a spectrum across the 1024 x 1024 pixel MAMA detector using the G140L 
grating for three different y-positions. For approximate wavelength correlations, the central wavelength of 
the grating is 1425 A (x=512), and the average dispersion is 0.6 A/pixel. 
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Fig. 8. — Spectral extractions utilizing our STISFIT and the standard CALSTIS routines for three stellar 
objects. N V, C IV, and several interstellar lines, marked with a *, are labeled in the middle panel. The 
STISFIT data are indicated by the heavy, black, solid lines, while the CALSTIS data are shown with lighter 
lines. 
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IMAGEfcy): 

Subtract background, B(x,y), from flat-fielded STIS image. 



cent(n): 

Locate center(s) of n stellar spectra to be extracted. 



FIT(x,Y,n): 
Build 2-D template of STIS image with n 
triple-Gaussian functions, correcting for variations 
in x and y using NORMJy), cent(n), and exttab. 



FINALFIT(x,y,n): 

Fit FIT(x,y,n) to IMAGE(x,y), 
modeling n spectra concurrently. 



SPECTRAfcn): 

Extract n net spectra from 
the grouped FINALFIT(x,y,n). 



2nd 



Run 



I 



FLUX(x,n): 
Flux-calibrate final spectra (ergs s _1 cm -2 A 1 ) 
by multiplying SPECTRA(x,n) by S A (x,n). j 



1 st Run 



shft(n): 

Calculate x (wavelength) offsets relative to center for 
each n using interstellar lines from SPECTRA(x,n). 



W/x,n): 

Create corrected wavelength arrays by using shft(n) 
to shift default wavelength array, W(x), for each n. 



S A (x,n): 

Create corrected sensitivity curves by mapping S(x) 
onto new wavelength scales, W A (x,n), for each n. 



Use shft(n) in 2 nd run of fitting routine to accurately 
align NORMj(y)W\th x-shifts when building FIT(x,y,n). 



Fig. 9. — Applying our developed models and software to the STIS dataset for NGC 604. (See Section 3.4 
for further detail.) 
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Fig. 10. — Comparison between our STISFIT low-dispersion extraction, CALSTIS low-dispersion extraction, 
and convolved CALSTIS high-dispersion extraction for NGC 346-368. The high-dispersion data are smoothed 
with a Gaussian to match the spectral resolution of the STISFIT data. The STISFIT spectrum is plotted with 
the thick line. The CALSTIS low-dispersion spectrum is plotted with the thin line and the high-dispersion 
spectrum is plotted with the dotted line. 
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Fig. 11 — STISFIT extracted spectra for the three closely spaced stars, 578A, 578B, and 578C, in NGC 604. 
(Also see Figure 2.) 
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Tabic 1. Hubble Space Telescope Datasets Used 



Object 


Instrument 


Filter/Grating 


Central A 

(A) 


Dataset Name 


GD153 


STIS 


G140L 


1425 


o4a502060 


NGC 346 


STIS 


G140L 


1425 


o46j01010 




STIS 


E140M 


1425 


o4wr01030 


NGC 604 


STIS 


G140L 


1425 


o4xl01040 




WFPC2 


F336W 


3344 


u2ab0207t 




WFPC2 


F336W 


3344 


u2ab0208t 




WFPC 


F439W 


4353 


wOnnOlOlt 




WFPC 


F439W 


4353 


w0nn0102t 




WFPC 


F439W 


4353 


w0nn0103t 


SAND 3 


STIS 


G140L 


1425 


o4xl01030 



